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ABSTRACT 

An accumulation of multi-epoch, high-resolution, optical spectra reveal that the nearby 
star G77-50 is a very cool DAZ white dwarf externally polluted by Mg, Fe, Al, Ca, 
and possibly Na, Cr, Mn. The metallic and hydrogen absorption features all exhibit 
multiple components consistent with Zeeman splitting in a i? « 120 kG magnetic field. 
Ultraviolet through infrared photometry combined with trigonometric parallaxes yield 
Toff = 5310 K, M — 0.60 M0, and a cooling age of 5.2 Gyr. The space velocity of the 
white dwarf suggests possible membership in the Galactic thick disk, consistent with an 
estimated total age of 8.6 Gyr. G77-50 is spectrally similar to G165-7 and LHS2534; 
these three cool white dwarfs comprise a small group exhibiting both metals and 
magnetism. 

The photospheric metals indicate accretion of rocky debris similar to that con- 
tained in asteroids, but the cooling age implies a remnant planetary system should be 
stable. A possibility for G77-50 and similarly old, polluted white dwarfs is a recent 
stellar encounter that dynamically rejuvenated the system from the outside-in. Metal 
abundance measurements for these cooler white dwarfs have the potential to distin- 
guish material originating in outer region planetesimals injected via fly-by. If common 
envelope evolution can generate magnetic fields in white dwarfs, then G77-50 and its 
classmates may have cannibalized an inner giant planet during prior evolution, with 
their metals originating in terrestrial bodies formed further out. Although speculative, 
this scenario can be ruled out if terrestrial planet formation is prohibited in systems 
where a giant planet has migrated to the inner region nominally engulfed during the 
post-main sequence. 

Key words: circumstellar matter — planetary systems — stars: abundances — stars: 
evolution — stars: magnetic fields — white dwarfs 



1 INTRODUCTION 

Metal-enriched white dwarfs have experienced a resurgence 
of interest coinciding roughly with the l aunch of Spitzer. 
Begin ning with the landmark study of iZuckerman et al^ 
(|2003l l. evidence has gradually accumulated supporting a 
picture whereby most, if not all cool, metal-lined white 
dwarfs obtain their photospheric heavy elements via rocky 
bodies similar in composition and in mass to large Solar 
System asteroids (|Farihi et al.ll2010al : |jurall2006l ). Data sub- 
stantiating this picture comes primarily from Spitzer photo- 
metric detections of infrared excess indicating circumstellar 
dust orbiting within the Roche limit of a substantial frac- 
tion of metal-contamin ated white dwarfs (|Farihi et al.ll2009l : 
Ivon Hippel et al]|2007l ). and via mid-infrared spectroscopy 
of the disk material, revealing silicate minerals ( primaril y 
olivines) associated with rocky planet formation (jJura et al.i 
l2009al : iReach et al.ll2005l ) . 



A powerful example of the scientific potential inher- 
ent in th ese stars is the spectacul a rly polluted whit e dwar f 
GD362 l|Kawka fc Venned l2005l : iGianninas et all 120041 '). 
This system shows remarkable infrared excess emission from 
closely orbiting dust, reprocessing over 3% of the incident 
stellar flux and exhibiting perhaps the strong est silicate 
emiss ion feature associated with any mature star (|jura et al.l 
l2007h . The optical spectrum of GD 362 displays 15 elements 
heavier than helium i n an abundance pattern mimicking the 
Earth- Moon system (jZuckerman et al]|2007^ . The convec- 
tion zone of the star has been enriched, at a bare minimum, 
by the equivalent of a 240 km asteroid, and possibly by a 
body as massive as Callis to or Mars, with promising ev- 
idence for internal water (|jura et al.l l2009"bl ). Thus, white 
dwarfs enable unique insights into terrestrial planet forma- 
tion around intermediate mass stars by providing data that 
can be obtained no other way; a lower mass limit to and 
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Figure 1. The normalized and coadded UVES spectrum of G77-50, rebinned onto a grid with spacing 0.1 A. The quasi-periodic pattern 
prominent at shorter wavelengths is due to a light path difference between the flat-field lamp and the sky; this interference is not fully 
correctable, especially at high S/N. Several previously unidentified metal absorption features are evident, including lines of Mg, Fe, Al, 
Ca, and perhaps Na, Cr, Mn. H/3 is detected for the first time in a very cool white dwarf. There are a large number of weak lines present 
in the spectrum whose authenticity is uncertain. The features in the green arm below 4700A and those in the red arm near 5900A and 
6250A are detector artifact and telluric absorption residuals (Table [TJ. 



the bulk chemical composition of extrasolar rocky, minor or 
possibly major planets. 

This paper investigates the spectral, kinematical, and 
atmospheric properties of the white dwarf WD 0322—019 
(G77-50). The star has been known as a very cool and metal- 
lined degenerate for over 35 years (|Hintzen fc Strittmattej 
1 19741 ). but was only recently found to exhibit split Can H 
and K absorptiori lines under high-resolution spectroscopy 
IZuckerman et all l2003l l. Despite similar, high-resolution 
spectroscopic observ a tions of several h undred white dwarfs 
l|Koester et al.l l2009l : IVoss et all 120071 '). including the de- 
tection of a few dozen with Call K absor ption features 
ijKoester et al.l l2005l : [Zuckerman et all l2003l l. only G77-50 
and possibly LTT8381 display evidence for Zeeman-split 
metallic lines. The observational program and scientific mo- 
tivation are described in §2, while the spectroscopic analysis, 
photospheric abundances, and stellar properties are derived 
in §3. The discussion of the results is contained in §4, where 
potential origins for the metals and magnetism in G77-50 
and related objects are explored in some detail. 



2 OBSERVATIONS & DATA 

2.1 Echelle Spectroscopy 

G77-50 was observed a total of 24 times between 2008 
October 4 and November 25 at Cerro Paranal with the 

8.2 m Very Large Telescope Unit 2 using the Ultraviolet and 



Visual Echelle Spectrograph fUVES: [Pekker et all \2md1 ). 

Spectroscopy was performed over the two detectors cov- 
ering wavelengths from 3200 A to 6650 A using a standard 
dichroic configuration with Ac — 3900/5640 A, resulting in 
two narrow gaps in spectral coverage near 4550 and 5650 A. 
A slit width of I'.'O was employed with 2x2 binning, re- 
sulting in a nominal resolving power of i? « 40 000 in both 
the UV-Blue and Red arms of the instrument. The obser- 
vations were taken at random intervals that were broadly 
logarithmic in temporal spacing, with each dataset consist- 
ing of two consecutive 900 s exposures. The featureless white 
dwarf WD 0000-345 (LHS 1008) was observed as a spectral 
standard on 2008 October 8 using an identical instrumental 
setup. 
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Table 1. Notable Features in the Coadded UVES Spectrum of 
G77-50 



Wavelength 


Absorbing 


Wavelength 


Absorbing 


(A) 


Element 


(A) 


Element 


3441 


Fel 


3934 


Can 


3570 


Fel 


3944 


All 


3582 


Fel 


3962 


All 


3619 


Fel 


3968 


Can 


3631 


Fel 


4227 


Cal 


3648 


Fe I 


'-tyj\j\j '-t 1 yjyj 




3720 


Fel 


4775 




3735 


Fel 


4861 


Hi 


3749 


Fel 


5073 




3758 


Fel 


5183 


Mgi 


3767 


Fel 


5890 


Nal 


3820 


Fel 


5900 




3832 


Mgi 


6250-6350 




3838 


Mgi 


6563 


Hi 


3860 


Fel 







Note. Owing to magnetism, the line centers cannot be determined 
accurately and the wavelengths given are approximate or labora- 
tory values. In some regions, many lines of Fel contribute to a 
single, complex feature but only one wavelength is listed for sim- 
plicity. Those features without a corresponding element are due 
(in part or in whole) to detector artifacts or telluric absorption 
residuals. 

The echelle data were processed with the UVES pipeline 
version 4.3.0, including cosmic ray masking, flat fielding, 
wavelength calibration, order merging, and distilled using 
optimal aperture extraction. Raw spectra produced in this 
manner had signal-to- noise (S/N) ratios that fell between 
16 and 22 at 5000A. The spectral standard data were inter- 
polated in wavelength to match the solution for the science 
target and smoothed by 121 pixels (3.5 to 4.0A). The science 
target was divided by the spectral standard and the result- 
ing spectrum was multiplied by an appropriate temperature 
blackbody to achieve a relative flux calibration. 

The 2008 UVES dataset was supplemented by two 
archival UVES observation s of G77 - 50 ta ken for the SPY 
survey in 2000 (Napiwot zki et al.l 120031 ). and fully re- 
duced Keck HIRES sp ectra obtained in 1999 and 1998 
llZuckerman et al ]|2003l) . After correcting for heliocentric ve- 
locity at the time of individual exposures, a master spec- 
trum was constructed from the normalized coaddition of all 
26 UVES spectra. The resulting spectrum has a bin width 
of O.lA, which increases the S/N, while reducing the origi- 
nal resolution to approximately match the decrease result- 
ing from the wavelength shifts observed in the Ca ll H and K 
lines (see §3). This spectrum is displayed in Figure [1] while 
Table [1] lists prominent stellar features and non-stellar arti- 
facts. 



3 SPECTROSCOPIC ANALYSIS & SYSTEM 
PARAMETERS 

G77-50 is listed as a binary suspect in lBergeron et akl (Il997l ) 
based on the shape of its Ha absorption feature and the lack 
of a trigonometric parallax distance at that time . Seeming 
to confirm this hypothesis. Zuckcrmar Tet al.l (|2003l ) reported 
the white dwarf as binary based on what appeared to be 



Table 2. Can H and K-line Velocity Measurements for G77-50 



HJD— 2 400 000 




v+ 


Instrument 


(d) 


(kms ) 


(kms 




51158.87475 


—43.79 ± 1.45 


-f82.10 ± 1.83 


HIRES 


51404.10902 


—51.36 ± 3.11 


+99.77 ± 3.61 


HIRES 


51737.93210 


-61.34 ± 5.43 


+94.93 ± 6.89 


UVES 


51740.88365 


—50.32 ± 4.75 


+87.94 ± 3.73 


UVES 


54743.71492 


—43.26 ± 2.69 


+98.32 ± 2.32 


UVES 


54743.72596 


—48.18 ± 2.71 


+92.99 ± 2.21 


UVES 


54743.73782 


—46.78 ± 2.78 


+96.59 ± 3.60 


UVES 


54743.74887 


—48.55 ± 2.61 


+91.46 ± 3.16 


UVES 


54743.81040 


—50.28 ± 2.87 


+91.03 ± 2.64 


UVES 


54743.82149 


—49.51 ± 2.15 


+86.41 ± 2.90 


UVES 


54744.73366 


—50.01 ± 1.79 


+89.36 ± 1.82 


UVES 


54744,74471 


—47.93 ± 2.43 


+91.32 ± 2.41 


UVES 


54744.83545 


—51.48 ± 1.99 


+89.90 ± 2.40 


UVES 


54744.84656 


—45.99 ± 2.04 


+87.58 ± 1.80 


UVES 


54747.69594 


—54.15 ± 2.32 


+97.44 ± 2.32 


UVES 


54747.70700 


—56.27 ± 2.51 


+94.08 ± 2.27 


UVES 


54749.65500 


—54.65 ± 4.22 


+88.54 ± 2.95 


UVES 


54749.66605 


—53.95 ± 2.97 


+93.02 ± 2.77 


UVES 


54750.66219 


-58.50 ±1.87 


+91.74 ± 1.96 


UVES 


54750.67321 


-52.67 ±2.09 


+92.08 ± 2.26 


UVES 


54759.63148 


-46.90 ±3.15 


+81.53 + 3.08 


UVES 


54759.64259 


-38.68 ±2.40 


+94.38 ± 2.40 


UVES 


54760.73176 


-44.87 ±2.32 


+79.17 + 2.55 


UVES 


54760.74281 


-42.05 ±3.26 


+85.12 + 3.15 


UVES 


54767.68013 


-44.27 ±2.07 


+76.15 + 2.15 


UVES 


54767.69118 


-38.95 ± 1.70 


+82.10 ± 1.35 


UVES 


54795.63988 


-41.59 ±2.30 


+81.30 ± 2.55 


UVES 


54795.65095 


-37.12 ± 1.92 


+79.29 ± 1.82 


UVES 



double lines of Ca ll H and K with large velocity separations, 
and possible changes in these velocities between two epochs 
of observation. On the basis of this sound interpretation, 
the original aim of the UVES program reported here was 
to obtain the orbital period and mass ratio of the putative 
binary by measuring and monitoring the velocities of the 
Ca II H and K line components. A typical, individual UVES 
spectrum of G77-50 is shown in the top panel of Figure [5] 
an d reveals the se emin gly d ouble Ca ll features first detected 
bv lZuckerman et ai.l (2OO3I ). 

3.1 Analysis and Interpretation of the 'Double' 
Calcium Lines 

3.1.1 Morphological Origins and Issues 

Initially, the two components seen in the Call H and K 
lines of each UVES spectrum were treated as originating 
in distinct stellar constituents in a spatially-unresolved bi- 
nary. New measurements taken in 2008 (MJD > 54000) were 
combined with velocities obtained via identical analysis of 
archival UVES and Keck HIRES data to determine an or- 
bital period. Each major component of the Ca 11 features 
were simultaneously fitted by a Gaussian plus a Lorentzian 
line p rofile with the programme FITSB2 (|Napiwotzki et aU 
I2OO4I ). The parameters of the line profiles were determined 
simultaneously for all spectra, i.e. assuming a constant line 
strength and shape (more on this below). The component 
velocities thus obtained are listed in Table [5] If the mea- 
sured variations in the line shifts were due to binarity, then 
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Figure 2. Upper: An example of an unbinned, single UVES spec- 
trum of G77-50 in the Ca II H and K region. Each individual expo- 
sure exhibits these double-peaked, metal absorption lines. Middle: 
The same region in the combined spectrum of all 26 UVES ex- 
posures. In this deep spectrum, broad absorption features due to 
All are revealed. Lower: In red and blue, the two coadded spec- 
tral subsets discussed in §3.1.1 are shown, plotted on top of the 
master spectrum in black. 



each stellar component should exhibit both red- and blue- 
shifted lines over time. Under this assumption, a sinusoidal 
period of almost exactly 2.0 d can be found within the line 
velocity data. 

However, if the observed changes in the Call H and 
K line components are the result of orbital motion in two 
components of a binary, then the resulting velocity ampli- 
tudes predict an essentially equal mass system. One can 




2,0 4.0 6,0 8,0 

Wavelength (jim) 

Figure 3. Spect ral energy distribution of G7 7-50 as revealed 
by GALEX llMart in c^al. 2005), UBVRIJH K llMcCook fc SionI 
ll999':'Berg eron et aljl997fl . and SpitzerlRAC ('Farihi et al.l2008f) 
photometry . The red error bars plot the observational data, and 
the black solid line is the stellar atmosphere model with param- 
eters listed in Table [3] 



then estimate the radii and masses of the components based 
on two inde p enden t trigonometric parallaxes for G77-50; 
ISmart et all (120031 ) report tt = 59.5 ± 3.2 mas or d = 
16.8to^ PC, while the US Naval Observatory (H. Harris 2010, 
private communication) measure tt — 58.02 ± 0.44 mas or 
d — pc. Taking the latter value by virtue of its 

smaller error, the predicted absolute magnitude of each 
equally luminous component would be My = 15.7 mag. At 
a nom inal effective temperature of 5200 K (|Bergeron et al.l 
Il997l ). this corresponds to 0.9 components. 

While this possibility is exciting as the total mass would 
then easily exceed the Chandrasekhar limit, it also leads to a 
major inconsistency. The ultraviolet through infrared pho- 
tometric energy distribution of G77-50 is consistent with 
a single effective temperature (Figure [3]), and thus if the 
system were binary, the stars would be twins to all limits 
of observation, in temperature, in mass-radius, and also in 
metal abundance as the Call H and K line components are 
equally strong. However, the profiles of the Ha line in in- 
dividual UVES spectra are found to be inconsistent with 
hydrogen-rich white dwarf models with log g > 8.1; high 
surface gravities fail to reproduce the observed line. The re- 
sults of the Ha line fitting are consistent with a single, cool 
white dwarf at d = 17.2 pc as determined by trigonometric 
parallax. 

Based on the discrepancies described above, the mas- 
ter spectrum was constructed as detailed in §2 and carefully 
examined by eye. The individual spectra exhibit only Ha, 
Ca II H and K prominently, but the combination of all UVES 
exposures yields a S/N « 100 spectrum with many addi- 
tional metal features and weak H/3 absorption. Importantly, 
all detected features in the coadded spectrum show evidence 
for multiple components, indicating the presence of a weak 
magnetic field. 

Given the spread in line component velocities as shown 
in Table [5J it is possible the lines will be smeared or other- 
wise diluted by the coaddition of all individual UVES spec- 
tra. To investigate the magnitude of this effect, the distri- 
bution of line velocities was searched for subgroups with 
similar values, whose coaddition might result in better line 
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sensitivity. Two groupings were identified and tiiese spec- 
tral subsets were coadded as above: 7 spectra with V+ in 
the range 75.5 — 84.5 kms~^ and 15 spectra with V+ in the 
range 87.5 — 95.5 km s~^. A section of these two spectra are 
shown in the lower panel of Figure [2] along with the master 
spectrum. The noise level of the subset spectra are increased 
compared to the master dataset, and there is no significant 
change in the shape or depth of the Ca ll H and K lines, nor 
the weak All lines. [It is noteworthy that the moderately 
strong Ca I line shown in the lower right panel of Figure |3] 
compares favorably in sharpness to both the non-magnetic 
and Zeeman-split models (§3.2).] 

3.1.2 Constraints on Line Variability and Rotation 

In this alternative picture, the velocity variation observed 
in the Call H and K line components is most likely due to 
stellar rotation. Assuming the star is not seen pole-on, the 
line of sight changes resulting from rotation span a (narrow) 
range of surface magnetic field strengths, and these induce 
changes in the strength of the observed Zeeman splitting. 

In order to search for periodicities in the wavelength 
shifts seen in the split H and K lines, a periodogram anal- 
ysis was performed for 1) AV = |V+ — Vl| and 2) V+ and 
V- independently. The first approach is adopted here, but 
the second yielded very similar results. A best fit is found 
for P = 29.84579 d and = 25.9, down from 280.5 when 
assuming constant line velocities, and thus indicating a high 
statistical significance of the variable solution. The phase- 
folded velocity differences confirm a good quality fit with 
Xred = 1.04, with a number of somewhat poorer fits are 
present in the period range 28 — 33 d. Two distinct solu- 
tions cannot be ruled out completely: 1.0311 d (x^ ~ 33.8) 
and 0.011333 d (x^ ~ 35.6). The latter period is suspect 
as it is similar to the length of the individual UVES ex- 
posures, while periods near 1.0 d are always dubious, and 
hence 29.85 d is the most likely solution. For a dipolar field 
configuration, half a rotation will produce a full cycle in the 
line velocities; if the long period solution is real, the white 
dwarf would be a slow rotator. 

The relatively strong Call features were examined for 
variation in strength over time within the UVES dataset. 
Measuring reliable equivalent widths for these blended lines 
was not possible, but the area swept out beneath the con- 
tinuum was examined as a proxy. Roughly half the feature 
strengths for both the H and K lines are found in a single, 
0.5(T-width peak, with all other measurements sitting below 
this peak, including some negative (i.e. unphysical) values. 
Based on this analysis, the individual spectra do not have 
sufficient S/N to meaningfully constrain variability in the 
observed line strengths. 

3.2 Stellar Parameters and Abundances 

When first studied. Ha was not detected in the optical spec- 
trum of G77-50 and the white dwarf was typed as DC in the 
classi fication scheme of the perio d (IHintzen fc Strittmatteij 
1974). Subseque ntly, ICreensteinl l| 19861 ) detected Ha and 
Sion et al] (|199G| ) typed the star in the modern spectral clas- 
sification as a DZA white dwarf; strongest lines are metallic, 
weaker line(s) of hydrogen. It is not clear whether feature 



Table 3. Stellar and Kinematical Parameters of G77-50 





5310 ± 100 K 


log g (cm s~^) 


8.05 ±0.01 


M 


O.eOiO.OlMg 


Mms 


1.4 ± 0.2 Mq 


^cool 


5.2±0.1Gyr 


tms 


3.4± l.SGyr 


ttot 


8.6± l.SGyr 


Ttrig 


58.02 ±0.44 mas 


(Ma > IJ-S ) 


(241, -877) masyr-i 


ZgC 


31.5 km s"'^ 




-11.2kms-i 




74.3kms-i 


V 


75.1 kms~l 


( U, V, W ) 


(-44.2, -60.5, -5.0) km s 



strength should be determined by depth or equivalent width 
(E. M. Sion 2010, private communication), but in the case 
of G77-50, a DAZ classification is readily arguable as its 
atmosphere is hydrogen-rich, with trace metals. 

As a first step, Tcft and log g are obtained by re-fitting 
existing BVRIJHK photometric data with the latest pure 
hydrogen model atmosphere grid and using the new paral- 
lax m easurement, following the method of iBergeron et akl 
(|l997l ). From this, an effective temperature of 5310 ± 100 K 
and a surface gravity of log g — 8.05 ±0.01 are ob- 
tained, corresponding to 0.60 ±0.01 Mq and a cooling age of 
5.2 ± 0.1 Gyr. Figure |3] plots the global energy distribution 
obtained with these stellar parameters (listed in Table ^ , 
together with ultraviolet through mid-infrared photometry 
obtained from both the ground and space. The consistency 
of the Spitzer data with the predicted stellar fiux indicates 
excellent agreement with the model, which was fitted only 
at shorter wavelengths. 

The coadded spectrum of G77-50 indicates the pres- 
ence of several elements as well as a magnetic field. From 
the splitting of the various atomic lines, B ~ 120 kG is es- 
timated. Such a field is too weak to have an important ef- 
fect on the thermodynamic structure of the star, and the 
absorption features are not strong enough to effect the de- 
termination of the atmospheric parameters described above. 
This was verified explicitly a posteriori by comparing a pure 
hydrogen model energy distribution with that of a blanketed 
model that takes into account the presence of all the heavy 
elements at the derived abundances. Hence, it is safe to use 
the above atmospheric parameters for the detailed abun- 
dance analysis presented below. It is, however, necessary to 
include the effect of Zeeman splitting in synthetic spectra 
calculations to obtain accurate heavy element abundances. 
As a first order approximation, magnetic line splitting is ex- 
plicitly included in the synthetic spectrum calculation by 
assuming the linear Zeeman effect in a constant magnetic 
field of 120 kC over the surface of the star. 

In the presence of a weak magnetic field, an atomic level 
with total angular momentum J splits into 2J±1 levels with 
magnetic quantum number m = — J, ...,±J. For each line, 
all transitions were calculated between the upper and lower 
levels that are allowed by the selection rule {rriu — mi = 0, 
±1, -1). Each transition is then shifted by 

AA = 4.67 X 10"'^ \l B (gimi - gumu) (1) 
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Figure 4. The combined UVES spectrum of G77-50 is shown in black and the best-fit model in red. Vertically ofli'set from the data, the 
same best-fit, but non-magnetic model is shown in blue for line identification. The following wavelengths correspond to their non-magnetic 
positions. Upper left: All Fel lines. Upper right: Pel lines and Mgl (3832,3838 A). Lower left: Call H & K, All (3944,3962 A), and Pel 
(near 3930 A). Lower right: Cal. 



where g is the Lande factor, B the magnetic field strength in 
G, m the magnetic quantum number of the level and Aq the 
central wavelength of the line in A. In the few cases where 
the Lande factor is not given in the VALeQ line list, the 
term designation is used to compute the Lande factor under 
the L-S coupling approximation, i.e. 



5 = 1 + 



JiJ + l) + SiS+l)-L{L + l) 
2J( J + 1) 



(2) 



The relative strengths of the vr, a (Am = 0, ±1 respectivel y) 
components are calculated according to ISobel'Manl (|l973l ). 

For the fixed value of T^a and log g determined above, 
a grid of synthetic, magnetic spectra is then calculated for 
each of the strongly-detected metals, as well as for more un- 
certain elements. The grids typically cover an abundance 
range from log [n(Z)/n(H)] = —8.0 to —11.0 in steps of 
0.5 dex. Next the abundance of each element is determined 
by fitting the various obser ved lin es using a similar method 
to that described in iDufour et al.i (j2005i") . The final adopted 

^ http://vald.astro.univie.ac.at 



abundance of an element is taken as an average of all the 
fitted features, and these are listed in Table [4] Finally, a syn- 
thetic spectrum including all the elements is calculated with 
the nominal abundances and plotted over the observed spec- 
trum in Figure [l] In order to facilitate the line identification 
of the various spectral features, a non-magnetic spectrum 
was also calculated with the same final metal abundances. 

Despite the simplicity of using a constant 120 kG field 
geometry for the synthetic spectrum calculations, the fits 
are remarkably good. Small discrepancies exist (e.g., in the 
wings and depth of the strong H & K lines), that are proba- 
bly an indication that the field strength varies slightly over 
the surface. It should be noted that the master, coadded 
spectrum represents an average taken at various, random ro- 
tation phases. Synthetic spectrum calculations using a more 
realistic (dipolar) magnetic field geometry, and a suitable 
average over the rotation period, would probably provide a 
better fit to the Call lines, but it is unlikely that such an 
effort would affect the derived abundances significantly. 

To find possible weak features arising from elements 
other than Mg, Fe, Al, and Ca, a synthetic spectrum was 
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Figure 4 - continued Upper left: Uncertain Mm features (4030,4033,4034 A), and Fel lines. Upper right: Uncertain Cri features 
(4254,4274 A), and Fel lines. Lower panels: Hydrogen Balmer lines a, 13. The dashed line in the Ha panel corresponds to a model 
100 K warmer than that used for the paper analysis (§3.3). 



calculated with all elements up to Ni in abundances s ome- 
what higher than their chondritic ratios (|Loddersll200^ ) rel- 
ative to the observed Ca in G77-50. This enhanced spectrum 
was laid over the coadded spectrum and a careful search was 
done by eye, allowing identification of the strongest possible 
features for additional elements as well as upper limits. The 
search revealed the possible presence of very weak features 
of Mm (4030,4033,4034 A) and Cri (4254,4274 A; see Fig- 
ure|4]). However, the putative features are only slightly above 
the noise level of the spectrum and thus these elements are 
considered uncertain. 

The Na I D doublet region also contains features (Figure 
[TJ, but caution is warranted as this portion of the spec- 
trum is contaminated by detector artifacts which signifi- 
cantly complicate the analysis. A careful examination of all 
26 individual UVES spectra indicates that the 5890 A line 
is real and thus Na is likely detected in this star. Unfor- 
tunately, the energy levels involved in these transitions are 
too close for the linear Zeeman approximation in a 120 kG 
magnetic field, and the Paschen-Back regime is appropri- 
ate. Hence, the position and strength of the lines cannot be 
matched with the calculations presented here, and an accu- 



rate analysis of Na in this star is best left for future work 
with cleaner data. Nevertheless, a rough Na abundance is ob- 
tained by matching the depth of the observed feature with 
the Zeeman approximation model, but interpretations based 
on this element should be avoided. 

Lastly, the synthetic Ha profile, appropriately split by 
a 120 kG magnetic field, is found to be slightly weaker than 
the observed line, whose strength depends sensitively on the 
effective temperature of the star. A model with a tempera- 
ture increased by 100 K produces a line profile that matches 
the observations well (Figure \^ , and falls within the effec- 
tive temperature uncertainty obtained from the fit to the 
photometric energy distribution. 



3.3 Stellar Kinematics 

Armed with an average, total velocity shift of 7 = 
-1-20.3 kms~^ from the median of all the Call line centers in 
Table [H a mass-radius constrained by trigonometric paral- 
lax, the stellar radial velocity of G77-50 can be disentangled 
from its gravitational redshift. The approximation 
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Table 4. Metal Abundances and Masses in G77-50 



Element ] 




(.10 g) 


tdiff 
(10'' yrj 


Nat 


-9.1 


0.25 


3.76 


Mg 


-8.3 


1.62 


3.85 


Al 


-9.2 


0.20 


3.72 


Ca 


-9.8 


0.08 


2.78 


Cr* 


-9.9 


0.08 


2.49 


Mn* 


-10.3 


0.03 


2.36 


Fe 


-8.7 


1.31 


2.33 


Total 


-8.0 


3.57 





(3) 



t Abundance uncertain, 
t Presence of this element is uncertain. 
Note. Errors are 0.2 dex 

GM 
^ c^R 

is valid for white dwarfs. For a star of 0.60 Mq and 
log g = 8.05, the radius is i? = 0.0121 Rq and Vg = 
31.5 kms^^. Thus the radial velocity of the binary system 
is Urad = 7 — Wg ~ — 11.2kms~^. Table 13] hsts the result- 
ing three dimensional space motion for G77-50 incorporat- 
ing the above radial v elocity and its observed proper motion 
l|Zacharias et al.ll2010l V While the velocities are not extreme, 
the V component lags behind the rotation of the Galaxy by 
60kms~^, indicating the white dwarf belon gs to the old thin 
or thick disk population (|Pauli et al.ll200^ ). 

A total age for the white dwarf can be estimated 
in the following way. Taking the average and standard 
deviatio n of three independent initial-to-fina l mass rela- 
tions ( WiUiams et a l. 2009; Kalirai et al. 2008; Dobbie et al.l 
I2OO6I ). the main-sequence progenitor of G77-50 had a likely 
mass of 1.4 ± 0.2 Mq. Main-sequence lifetimes for stars in 
this mass range ar e obta ined from the analytical formu- 
lae of iHurlev et al.l (|2000l V and added to the cooling age 
of the white dwarf. This results in a total system age of 
8.6 ± 1.3 Gyr for G77-5Q, a range consistent with thick disk 
membership, although not conclusively. 



4 ORIGIN OF METALS AND MAGNETISM 

Some of the material presented in this section is necessarily 
speculative, and are first attempts to understand how a star 
such as G77-50 came to possess its somewhat unusual char- 
acteristics. It should be the case that near-future theoretical 
and empirical investigations will be able to test these ideas 
more rigorously. 



4.1 Nature of the Photospheric Metals 

Metal diffusion timescales and convection zone parameters 
for G77-50 were calculated using a complete model atmo- 
sphere with the parameters and heavy element abundances 
from the preceding analysis in the outer layers. This was 
continue d to the bottom of the convection zone using pure 
hvdrogen lKoesteil (|2009h . Table[4]lists the relevant lifetimes 
for all the metals detected or suspected in G77-50. The 
total mass of the mixing layers in G77-50 is found to be 
1.23 X 10^^ g, and can be combined with the photospheric 



abundances to obtain the current mass of each heavy ele- 
ment within the star. As the metals continuously sink, these 
calculations yield the minimum mass of each accreted metal, 
and are listed in Table |3] The total mass of metals residing 
in the star is 3.57 x 10^" g and equivalent to that contained 
in a 65 km diameter body of density 2.5 gcm"'^; the highest 
mass of heavy elements yet inferred within a hydrogen-rich 
white dwarf. 

Owing to its very cool effective temperature, G77-50 
possesses a sizable convection zone on par with t hose in- 
trinsic to helium-rich white dwarfs (iKoestedlioogl '). Thus, 
the star has a commensurately long diffusion timescale 
for heavy elements; e.g. 280,000 yr for Ca. In this and 
other respects, G77-50 is characteristic of the coolest metal- 
rich white dwarfs. It has modest metal abundances, with 
log [n(Z)/n(H)] < -8.5 for all but Mg, a relatively low 
time-averaged accretion rate of 3.7 x 10^ gs~^ for all clearly 
detected and suspected heavy elements, and a lack of in- 
frared excess as measur ed with Spitzer IRAC photometry 
dFarihi et al.ll2009l . [20oi ). indicating a dearth or absence of 
dust in its immediate circumstellar environment. 

The relatively old white dwarf has a total space ve- 
locity of 75kms~^ and its current location is 17.2 pc 
from the Sun. At this speed, the star travels nearly 
77pcMyr~^ and may be passing throug h the relatively in- 
terstellar matter-deficie nt Local Bubble l|Redfield fc Linskvl 
l2008l : IWelsh et al.lll999l ), its metal content decaying from a 
past event within a dense region of gas and dust. Within a 
period just over 1 Myr, G77-50 could have been outside the 
Local Bubble, and its metal abundances would have decayed 
by a factor between 20 and 100, depending on the element. 
With roughly 2 x lO'^'^ g of atmospheric metals in this pic- 
ture, how dense an interstellar region is necessary to account 
for the accretion in such a manner? For interstellar heavy 
elements typically contained in dust p articles, the accretion 
rate onto a star is (|Farihi et al.ll2010al ): 



Mz = 



TiGMRp^ 



Tcv 



(4) 



where M and R are the stellar mass and radius, s — 
y/v'^ + Cs^ is the space velocity v in the supersonic regime 
(v ^ Cg) appropriate for G77-50, and Tev is the evaporation 
temperature of the dust grains. Conservatively estimating 
that dust evaporates at 1000 K, the minimum density re- 
quired to deposit 2 x 10^^ g of metal in the star over a typi- 
cal sinking timescale of 3 x 10^ yr is 1600 cm~^. While such 
densities exist in molecular clouds, the space velocity of G77- 
50 requires continuous accretion over a distance of 30 pc in 
order that sufficient material is captured. This scenario is 
unlikely; one of the nearest and largest cloud complexes, the 
Orion Nebula is only about 8 pc in diameter. 

Alternatively, the metals in the star were accreted from 
its circumstellar environment and were originally contained 
within remnant planetary bodies rich in heavy elements. The 
detection of both Al and Ca in the star supports a scenario 
in which the metals originated in a refractory-rich source 
such as a rocky, minor or major planet. This possibility is 
likely based on the weight of evidence favoring circumstel- 
lar accretion among the observed p opulation of single , cool , 
and m ctal-cnriched white dwarfs jZuckerman et al] I2OI0I : 
iFarih i ct al. 2010a). In this scenario, the photospheric met- 
als in G77-50 were delivered by one or more rocky planetary 
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Table 5. Metal Ratios in G77-50 and Other Polluted White 
Dwarfs with Al 



Star 


Name 


Mg/Fe 


Al/Fc 


Ca/Fe Disk? 


0322-019 


G77-50 


2.9 


0.31 


0.08 






1.2t 


0.14t 


0.06t 




solar 


1.2 


0.10 


0.07 


0208+096 


G74-7 


1.0 


0.08 


0.11 


0300-013 


CD 40 


1.7 


0.21 


0.40 + 


1633+433 


G180-63 


2.1 


0.14 


0.18 


1729+371 


CD 362 


0.5 


0.18 


0.26 + 



^At the end of an assumed accretion event 0.5 Myr prior. 

Note. The three middle columns list the relative number abun- 
dances. Metal ratios for oth er white dwarfs with Al detections 
are taken from the literature l lKlein et al.ll2010l : [ Zuckerman et al.l 
|2007|. l2003l : lLodderj[20o3'l . 

bodies passing close enough to the whi te dwarf to become 
tidally destroyed or otherwise accreted (|jurall200^ ). 

For a prior pollution event, the current metal-to-metal 
ratios have been altered from the accreted values, via the in- 
dividual heavy element sinking timescales. Table [5] lists the 
photospheric abundance ratios, relative to Fe, for the four 
heavy elements confidently identified in G77-50. Notably, 
both Mg/Fe and Al/Fe appear enhanced relative to solar 
values and to the stars with circumstellar dust; precisely as 
expected in a declining phase because Fe sinks most rapidly. 
Interestingly, Ca/Fe is nearly solar, yet down by a factor of a 
few compared with the stars currently accreting from disks. 
Under the questionable assumption that the accreted abun- 
dances were close to solar, then the photospheric pollution 
halted roughly 0.5 Myr ago. In all likelihood, the accretion 
epoch for G77-50 ended within the last few to several diffu- 
sion timescales, no longer than a few Myr ago. Events older 
than 3 Myr ago imply an Fe-dominated (Mpc/Mz ^ 94%) 
parent body more massive than Pluto. 



4.2 A Late Instability Near 5 Gyr 

Again, G77-50 epitomizes the older and cooler, Tofi < 9000 K 
metal-rich white dwarfs: something weighty has occurred re- 
cently in these Gyr old, presumably stable systems. Gener- 
ally, the timescale for a planetary system to dynamically 
settle - pre- or post-main sequence - should be 100 Myr 
l|Debes fc Sigurdsson|[20o3 ). and thus a catastrophic event 
at Gyr epochs is not expected for G77-5Q and similar white 
dwarfs. Dual planet interactions will occur during this short 
period, if they occur at all, but instabilities among three (or 
more) planets may occur on longer timescales, as hypoth- 
esized for the period of Late He avy Bombardment in the 
Solar System (iGomes et al-lbOOSl). 

Following Debes fc SigurdssonI ||2002D and 
19961 ) for the simple case of three 



[chambers et al.l | 



l|2002h 

case of 



0.001 Mq planets, one can calculate system architectures 
that give rise to an instability timescale of 5.2 Gyr during 
the white dwarf evolutionary stage. Taking initial and 
final stellar masses to be 1.8 and 0.6 Af0 produces an 
(adiabatic) orbital expansion factor of 3 between the main 
and post-main sequence phases. For an innermost planet 
now located at 10 AU, the outer planets would be near 30 



and 157 AU, with original - and pres umably m etastable - 
semimajor axes of 3.3, 10, and 52 AU jPebes fc SigurdssonI 
|2002| ). In general, for all layouts that keep the innermost 
planet safely outside the stellar photosphere (ai > 2AU) 
during the asymptotic giant branch, the planets must be 
spaced such that a^jax > 3 and azja^ > 10. Thus, the 
onset of instability at Gyr timescales requires a somewhat 
finely-tuned, increasingly wide planet spacing. 

An alternate way in which an old and stable planetary 
system might be efficiently agitated is a stellar encounter. 
Though rare, at the cooling age and space velocity of G77- 
50, it can be shown that at least one fly-by is consistent with 
theoretical expectations. Dynamically, the number of stellar 
encounter s per unit time, within a d istance D of a star is 
given by (|Garci'a-Sanchez et al.lll999l ): 



N — TvD v,pt 



(5) 



where Vt is the space velocity of the white dwarf relative to 
passing stars, and p, is the local density of stars and stel- 
lar systems. From Hipparcos data, the velocity dispersion of 
stars within 2 kpc of the Sun and |fe| > 30 ° v aries between 2 2 
and 44kms~^ depending on spectral type (|Mignardl |200G| ) . 
Later type stars tend to be kinematically more stirred as 
they represent, on the whole, relatively older populations. 
And because the most abundant stars in the Galaxy are K 
and M dwarfs, these relatively fast moving stars are the ob- 
jects for which encounters are most likely. Taking 40kms~^ 
for a typical field star and the 75kms~^ space velocity of 
G77-50 yields Vt = 85kms~^. The local space density of 
stars is 0.081 pc~ ^ (T. Henry 2010, private communication; 
iHenrv et all2006l fl and Equation [5] gives N = 22 Myr"^ for 
encounters within 1 pc. For close fly-bys that might signifi- 
cantly impact the outer regions of a planetary system, take 
D ^ 1000 AU and then = 0.5Gyr"\ If correct, and such 
an encounter is 50% probable within each 0.5 Gyr window, 
then it is 99.9% likely that G77-50 has suffered at least one 
such stellar encounter within its 5.2 Gyr cooling history. 

An encounter could disturb a Kuiper Belt analog origi- 
nally orbiting at a few to several tens of AU, but now resid- 
ing at one to a few hundred AU owing to mass lost during 
the post-main sequence. Large objects such as Sedna, whose 
mass is similar to Ceres, can readily be perturbed into high 
eccentricity orbits by a stellar encounter , thus sending outer 
plane tesimals towards the inner system (jKenvon fc BromlevI 
I2OO4I ). While chances are tiny that any single, distant body 
could be flung directly within the tidal disruption radius of 
the white dwarf, the perturbation of many bodies (i.e. a sig- 
nificant fraction of a surviving population) could result in 
their capture and further s cattering within the inner system 
(|Debes fc Sigurdsson|[20o3 ). 

A potential hurdle for this hypothesis is the total avail- 
able mass that survives heating during the giant phases at 
these orbital distances. Out to roughly 100 AU, objects com- 
posed of pure water ic e and up to 100 km in diamet er should 
sublimate completely l| Jurall2004 IStern et af]|l99Gl '). but the 
overall effect on realistic bodies of heterogenous chemical 
composition is unknown. For example, some models predict 
that subsurface volatiles should be protected by superior 



^ Statistics of the solar neighborhood are continually updated at 
[http : / / www . recons . org | 
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layer s of non-volatile material such as silicates (|Jura fc 
l201Cf ). Studies of Kuiper Belt analogs in the post-main se- 
quence may reveal the overall impact of sublimation and 
constr ain the total mass tha t survives to the white dwarf 
stage (|Bonsor fc Wvattll2010l ). 

In contrast, Oort cloud analogs will not be destroyed 
via heating but are prone to dynamical evaporation in the 

! iost-main sequence if th e stellar mass loss is asymmetric 
Parriott fc AlcockllT998l '). At tens of thousands of AU, stel- 
lar encounters can be frequent over Gyr timescales typi- 
cal of very cool white dwarfs such as G77-50. These events 
should strip away some portion of any Oort-like cloud, but 
also perturb another fraction onto eccentric orbits overlap- 
ping with the inner system. If such cold, outer planetesimal 
belts are ultimately responsible for some fraction of metal- 
polluted white dwarfs, then a large deposition of volatile 
elements would be expected, but h as not yet been observed 
l|Zuckerman et al.ll2007l : |jurall2006l ). 

The subset of older and cooler metal-rich white dwarfs 
that G77-50 represents are a relative ly high velocity group o f 
stars, typically with T > 50kms~^ (jAannestad et al.lll993h . 
and the above exercise can be broadly applied to Gyr age 
white dwarfs with metals. Thus, stellar encounters have the 
potential to make ancient planetary systems dynamically 
young for a brief period, and may account for the popula- 
tion of very cool DAZ and DZ stars, of which vMa2 is the 
prototype. 

4.3 Weak Magnetism 

In the favored model for cool, metal-rich white dwarfs, a 

J ilanetary system has surviv ed post-main sequence evolution 
Debes fc SigurdssonllioO^ V Instabilities drive rocky plane- 
tary bodies such as asteroids into close approach with the 
compact star, where they become tidally- or otherwis e de- 
stroy ed and are sometimes observed as dust disks (|jural 
|2003| ). For the high eccentricity orbits necessary for the 
(eventual) delivery of metals onto the surface of the star, 
conventional sized planets are the most efficient perturbers. 
In this picture, an otherwise-replete planetary system per- 
sists at metal-enriched white dwarfs, likely truncated near 
the maximum radius of the asymptotic giant progenitor star. 

Based on precision radial velocity planet searches, 
10.5% of FGK dwarfs with masses in the range 0.7— 1.3 M p, 
host gas giant planets within 3AU jCumming et al.ll2008h . 
Due to the finite time period over which these searches have 
been carried out, most of these giant planets orbit within 
1 AU, and a decent fraction orbit close to their host star 
(the so-called hot jupiters). Relative to these findings for 
solar-type stars, giant planets are found more frequently, and 
with a higher mass distribution, in radial veloci t y sea rches 
of intermediate mass stars. Both lLovis fc MavoJ (20o3) and 
Ijohnson et al.l |20o3) find that planet frequency and planet 
(minimum) mass scale with stellar host mass, and that the 
fraction of gas giants orbiting within 3AU of stars with 
masses between ar ound 1.5 and 2.0 M© is closer to 25% 
ijBowler et al.ll2010l ). Interestingly, stars with higher masses 
do not host closely-orbiting giant planets, possibly due to 
halted orbita l migration from rapid inner disk dis sipation 
JCurriellioogl '). or planet-star gravitational tides jHansenl 
I2OIOI ). 

The bulk of giant planets orbiting within a few AU of 



their host stars will be destroyed during the post-main se- 
quence phases of their host stars. During the first ascent gi- 
ant branch (RGB), stars with main-sequence masses above 
1.0 Mq expand to around 10 - 20 ii© , or 0.05-0.10 AU, suf- 
ficient to swallow only the hot planets (| Bowler et al.|[2O10l : 
IVillaver fc Livid l2009i v Depending on which initial-to-final 
mass relation one uses, the main-sequence progenitor of G77- 
50 was an early to late F star with a mass between 1.2 
and 1.6 M0. For a M > 1.3 Mq progenitor, there is a 25% 
chance it hosted a giant planet within a few AU, but not 
a hot planet. In this case, the inner planet would survive 
the RGB intact but most likely become enveloped during 
the asymptotic giant branch (AGB), either directly or via 
tidal for ces; a 1.5 Mp, star has a maximum AGB radius near 
2.4 AU (|Villaver fc Lividl2007h . 

Depending on the mass and angular momentum of an 
engulfed planet, a common envelope may d evelop for a time 
suffic i ent to generate a magnetic dynamo ( Nordhaus et al. 
2010*: 'To ut et al.ll2008l : ISiess fc Liviolll999bl : iRegos fc Tout 



1995). The eventual accretion of the planet onto the stellar 
core may also result in enhanced mass loss from the giant 
(|Siess fc Livid[l999a^ . producing a white dwarf remnant po- 
tentially less massive than predicted by initial-to-final mass 
relations. If this occurs at the level of 0.05 Mq in the white 
dwarf, then the inferred mass of the progenitor star will be 
bi ased towards lower, m ain-sequence masses. In the model 
of lPotter fc Tb^ (|2010l ). the longer the common envelope 
environment, the larger the magnetic field strength in the 
resulting white dwarf plus companion merger. This model 
predicts fields up to 10^ G in a common envelope lasting 
roughly 10"* yr, but with a stel lar companion genera ting the 
convective motions. Following IPotter fc ToutI (|201G| ) but re- 
ducing the companion mass by a factor of 100 (to 10 Mj), 
all else being equal, a magnetic field of 10^ G is predicted. 
Interestingly, the surface field of G77-50 is 1.2 x 10^ G. 



4.4 Possible Planetary System Architecture 

In the admittedly speculative picture outlined above, the 
primordial G77-50 system formed a gas giant that mi- 
grated through a disk to an inner orbit. Subsequently and 
bey ond this orbit, additional planets may have formed 
(.Mandell et al.ll2007l ). As in the Sola r System and a signifi- 
cant fraction o f main-sequence stars (|Carpenter et al ] |2009l : 
ISu et al.ll2006l ), a Kuiper belt analog existed at G77-50, con- 
sisting of primordial icy and rocky bodies condensed from 
the stellar nebula. During post-main sequence evolution, 
the volatile outer layers of these planetesimals were lost by 
evaporation and radiation pressure, leavi ng behind objects 
with non-volatile surfa ces and intact cores l|Bonsor fc WvattI 
I2OIOI : lJura fc Xull2010^ . 

In order that a stellar encounter repopulate the in- 
ner, circumstellar regions around any star, sizable planets 
must exist there to trap or further scatter injected bod- 
ies. Such captured planetesimals are vulnerable to perturba- 
tions via newly establ i shed r esonances on 100 Myr timescales 
(|Debes fc Sigurdssod I2OO2I ). From this point forward in 
time, the eventual pollution of the star by a tidally destroyed 
planetesimal proceeds as envisioned for metal-rich white 
dwarfs with cooling ages less than 0.5 Gyr (jParihi et al.ll2009l : 
[jura 2003). If the stellar encounter hypothesis applies to 
G77-50 and other very cool, polluted white dwarfs, then the 
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atmospheric metals may be a reflection of the rocky remains 
of primitive, outer system planetesimals rather than Solar 
System asteroid analogs. As such, they would have a com- 
position distinct from objects that formed at higher temper- 
atures in the inner system, presumably including a higher 
concentration of volatile elements. 

Testing this hypothesis may be difficult. While cooler 
white dwarfs have lower atmospheric opacities and hence de- 
tection of their heavy elements is less challenging, there are 
two complications. The main difficulty is that the coolest 
DAZ and DZ stars are not amenable to ultraviolet spec- 
troscopy, which is the ideal region to search for volatiles 
(or any element) in low abundance. For example, trace 
amounts of atmospheric carbon are only seen in the opti- 
cal spectra of cool, helium-rich whi te dwarfs, and only at 
abundances log [n(C)/n(He)] > -7 (|Dufour et al.ll2005l ). If 
C/Fe — 8.3 in G77-50 as in chondrites, then it would have 
log [n(C)/n(H[)] = —7.8 and be undetectable in this manner. 
Also, it is not clear what fraction of volatile material would 
remain after heatin g and ablation during the post-main se- 
quence (jura 2006") . 

An additional complication is that it is difficult to con- 
strain the epoch of accretion in these cooler, metal-polluted 
stars. Observationally, it is rare that cool DAZ or DZ have 
dust disks, and this may be the result of a significant reduc- 
tion in the ratio of disk lifetime to cooling age. The observed 
abundances can be used to constrain the epoch of accretion, 
but this is fundamentally biased towards the composition of 
Solar System objects and may not accurately reflect those 
of extrasolar planetary bodies. This latter concern may be 
partially alleviated in the growing num ber of DZ stars bein g 
found in the Sloan Digital Sky Survey (|Dufour et al.|[2007i ): 
follow up of these stars will better constrain the relative 
frequency of disks at older, metal-enriched stars and their 
abundances will hopefully provide clues to the nature of the 
polluting bodies. 



5 CONCLUSIONS 

The cool white dwarf G77-50 is found to be simultane- 
ously metal-contaminated and magnetic. A hypothesis that 
may account for both properties is an evolved planetary 
system. While high ly uncertain, it is possible that G 1 65-7 
(Dufour et all bood l and LHS 2534 (|Reid et al.l [2OOII ) ac- 
quired their magnetic flelds and photospheric metals in 
a manner analogous to G77-50. In this tentative picture, 
an inner giant planet was swallowed by the post-main se- 
quence progenitor, creating a common envelope whose con- 
vective motions generated the extant magnetic field. Sur- 
viving rocky planetesimals have polluted the star with sev- 
eral refractory metals that are typical major cons tituents 
of as teroids and inner Solar System planets ( AUcgre et al.l 
I1995I ). However, with a cooling age of 5.2 Gyr, it is difficult 
to imagine a planetary system architecture with instabili- 
ties sustained over such a period. A plausible model is a 
stellar encounter that perturbs outer planetesimals into the 
inner system, where planets scatter and trap them with a 
renewed dynamical potential for close approaches with the 
white dwarf. 

The exact nature of the parent body whose metals 
currently pollute G77-50 will likely remain uncertain for 



the foreseeable future. The white dwarf has been observed 
over a large wavelength range with the two most power- 
ful, high-resolution optical spectrographs currently avail- 
able, and with the exception of Na, the detections and abun- 
dances here are not likely to be improved upon. Further- 
more, the effective temperature of the star means that there 
is little flux to observe with the ultraviolet spectrographs on 
HST. G77-50 and other very cool DAZ and DZ stars will 
have to await more powerful long wavelength observations 
with JWST or ALMA. 
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